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Solar flares presumably have an impact on the deepest layers of the solar 

^— ^ I atmosphere and yet the observational evidence for such an impact is scarce. 

Using ten years of measurements of the Na Di and Na D2 Praunhofer lines, 

j^i measured by GOLP onboard SOHO, we show that this photospheric line is 

indeed affected by flares. The effect of individual flares is hidden by solar 

oscillations, but a statistical analysis based on conditional averaging reveals 

^~: ' a clear signature. Although GOLP can only probe one single wavelength at a 

time, we show that both wings of the Na line can nevertheless be compared. 

■^ ■ _ The varying line asymmetry can be interpreted as an upward plasma motion 

p^ ■ from the lower solar atmosphere during the peak of the flare, followed by a 

• . downward motion. 

O ' 1 Introduction 

u. 

52 ■ Solar flares are the most powerful events occurring in the solar system. The mechanism 

invoked for their energy release is the conversion of magnetic energy into radiation, ther- 
mal and kinetic energy through the reconnection of magnetic field lines. More precisely, 
^ ' the coronal magnetic field suddenly changes its configuration and free energy when the 

£r) . stre ss imposed at its footpo ints by photospheric plasma motion exceeds a certain thresh- 

\^ [ old (JArchontis et all l2005l). The deta ils of why and how the energy is released are still 



o 



en ' largely debated {e.g ( Schriived . l2009l )). Once converted, the kinetic energy accelerates 

l/^ ■ particles; part of these are directed downwards and deposit their energy in the solar 

^P ! atmosphere while another part can eventually escape into interplanetary space. The 



interaction of these accelerated particles with the surrounding plasma ultimately results 



in ele ctromagnetic radiation at various wavelengths, from X-ray to radio domains (|Bena . 
200i). 



rS I One of the best observed manifestations of flares is the appearance, after the im- 

C^ ' pulsive phase, of loops that emit strongly in the extreme ultraviolet and soft X-rays 

(SXR). Downward-accelerated particles lose their energy by collisions when they en- 
counter the dense chromospheric plasma. The local material is heated, expands, and 
then rises to fill the newly configured magnetic-field loop where plasma cools down 
thro ugh radiative and conductive los s es. This is the mod el of chromospheric evapora- 
tion ( Antonucci. Gabriel, and Dennisl . ll984l : IPisheii . Il987bl ): the hot plasma located just 



above the energy-deposition layer evaporates towards the corona. The cool plasma just 
underneath is pushed down; this is called chroniospheric condensation. In this model, the 
energy deposition layer is located in the middle chromosphere or in the transition region. 
The rise of chromospheric plasma is manisfested by the observation of Doppler-shifted 



spect ral lines that form at these temperatures ([Schmieder et all 119871: iBerlicki et al. 



20051 ) and by the downflow of plasma as observed from the Ha line (JFisheiJ . Il987al ) . 



White-light (WL) continuum emission has been also observe d in some flares - that are 

then called white light flares ( WLF) - but their origin is unclear (JHudson. Wolfson. and Metcali 
20061 ). It is generally agreed that the WL emission takes place in the lower atmo- 



sphere, near the minimum temperature region or below ( Aboudarham and Henouxl . 
19861 ; iDing et al. , 2003). How these layers are heated, however, remains unknown (see 
the review by (JDingi . 120071 )). Among the various proposed mechanisms, direct heating by 
non-thern ial electron beams, chromospheric radiat i ve ba ck-warming, or H~ continuum 
emission ( Xu et al\ . l2006l : iGan. Henoux. and Fand . I2OO0I ) are often cited. When taken 
individually, none of these mechanisms can provide the amount of energy required by 
the observa- t ions. Direct he a ting i n or near the photosphere has also been suggested 
( Ding et all Il994l : ISchrijveii . l2009l ). The latter mechanisms suggest that the low so- 
lar atmosphere (the chromosphere of course, but also the temperature minimun region 
and the upper photosphere below) is heated during these flares. As technology evolves , 
more WL emission is observed during flares ( Jess et aU I2OO8I : iKretzschmar et aUl2010l ). 
We present here observations in visible light made by the Global Oscillations at Low 
Frequency (GOLF: Gabriel et al, 1995) instrument onboard SO HO (space mi s sion o f 
international cooperation between ESA and NASA) since 1996 (JGabriel et all 119971 ). 
These observations give information on the low solar atmosphere during solar flares. 



One question that we address in this paper is the recurrent problem of retrieving a 
global picture of the flare since, in addition to the intrinsic variability of the events, the 
conditions (position, amplitude, ...) vary and the observations are often made by differ- 
ent instruments. We follow here the same approach as in Kretzschmar et al, (2010) and 
search for a statistical signature of flares by analysing ten years of SOHO/GOLF data 
consisting in the intensity of the Fraunhofer absorption Na D lines integrated over the full 
solar disk. T hese lines are formed in the so l ar photosphere below the temperatu r e min - 
imum region ( Caccin. Gomez, and Robertil . Il98d : iBruls. Rutten. and Shchukinal . Il992l ). 
We refer here to the photosphere as a few hundred kilometer thick region between the 
altitude where the optical depth at 500 nm is unity (T500 = 1) and the temperature 
minimum region. More precisely, the wings of the Na D line that we study in this paper 
are mostly gen erated between 30 and 500 km, just below the transition region and the 
chromosphere ( Elbe et al\ . I2OOII ). Simulations and observations for a quiet Sun atmo- 



spher e have been recently made for the formation of the Na D lines ([Leenaarts et al. 



2OIOI ): this study strongly contributes to the idea that these lines are originating from 
the photosphere. However during a flare, the formation heights of the Na D sodium may 
differ from quiet conditions. Nevertheless changes in the Na D lines during flares thus 
indicate that flares affect the lower solar atmosphere including both photosphere and 



chromosphere. Our analysis therefore provides an opportunity to study the influence of 
solar flares on the lower solar atmosphere. This is particularly interesting since the best 
observed flare features usually occur in the upper solar atmosphere, where the contrast 
between the flare signal and the background emission is largest. 

The intensity of the Sodium line has been measured every five seconds by GOLF. 
The good cadence, duration and duty cycle of the observations, together with the high 
quality associated with space-based observations, partly allow us to compensate for the 
lack of spatial resolution. As we shall see below, large and even medium amplitude 
flares do have a significant impact on the Na line, which we interpret as the signature 
of photospheric heating. Furthermore, this work also points to the possible influence of 
flares on helioseismic studies. 

In Section 2 we introduce the GOLF instrument, the data, and the analysis technique. 
In Section 3, the results are presented and discussed. Conclusions follow in Section 4. 

2 Instrumentation and Data Analysis 

2.1 The GOLF Instrument 

GOLF, which is one of the three helioseismic instruments onboard SOHO, is a resonant- 
scattering spectrophotometer that monitors periodic variations of the solar surface by 
observing the Doppler shift of the Na Di and Na D2 Fraunhofer lines respectively at 589.6 
and 589 nm integrated over the solar disk. The GOLF instrument cannot distinguish 
these two lines and so it monitors the sum of the signals coming from both lines. These 
two lines have approximately the same behaviour, and so we can have an idea of the 
evolution of one absorption line even when working with the sum. The principle of this 
device, however, is best understood by considering one single line. 

GOLF was designed to measure the solar flux in the wings of the absorption line, at 
± 0.0108 nm from the central rest wavelength as shown on Figure [TJ By comparing the 
intensity of the two wings, the Doppler shift can be retrieved and the frequency of solar 
surface oscillations can be determined. 

GOLF initially measured alternatively the intensity of each wing every 20 seconds. 
Unfortunately, due to a malfunction of the device, a single-wing working cycle had to be 
adopted to ensure the safety of the instrument. We consider here red-wing observations 
from 15th, October 1998 to 17th, November 2002 and blue-wing observations from 11th, 
April 1996 to 28th, June 1998 and from 18th, November 2002 to 23rd, May 2006. Note 
that an observed blue shift (negative velocity by convention) corresponds to a motion 
towards the observer, i.e. upward in the solar atmosphere. In the following, the SXR 
flux is measured by the SEM instrument onboard GOES in the 0.1 - 0.8 nm band. 




-20 -10 10 

Doppler Velocity (Km/s) 



Figure 1: Location of the two wavelengths measured by GOLF. When working with 
Doppler shifts, one commonly expresses wavelengths in terms of Doppler ve- 
locities. The measured intensity corresponds to it 5.5 km s~^ around the 
wavelength at rest. 



2.2 The Method: Conditional Averaging 

The signature of a single flare has no discernible effect on the intensity measured by 
GOLF because it is masked by omnipresent temporal fluctuations at the photospheric 
level, which are caused by internal acoustic waves and by granulation. GOLF was 
designed to study these pressure waves, but for us they are background noise that mask 
the flare signal. To overcome this problem, we refrain from looking at individual events 
and consider the flare signature in a statistical sense. We use a conditional-averaging 
technique in order to see the effects of solar flares. Th is technique has fr equently been 
used to detect coherent structures in turbulent flows ( Block et all 120061 ) and is called 
superposed epoch analysis in the frame of geophysics. By averaging N time sequences 
that are synchronised with respect to the same trigger (in our case the occurrence of 
the flare, as evidenced by the peak intensity of SXR flux), incoherent fluctuations are 
averaged out whereas reproducible features such as flare signatures add up. 

We sorted the flares in decreasing order of amplitude, from the most powerful X-class 
flares down to the weaker M, C, and B-class flares, by using the soft X-ray classes derived 
from GOES observations. The largest flares are often accompanied by solar energetic- 
particle events that perturb the instruments. To avoid a possible impact on our results, 
we discarded flares whose peak energy exceeded 5 x 10~^ Watts m~^ and checked that 
the remaining flares were not accompanied by such particle events. We excluded in that 
way the ten largest flares out of a total number of 116 X-class flares. Table [1] indicates 
the number of events observed in each wing. 



Table 1: Number of events for different class-flare and the wing in which it is observed 



X M C 



B 



Red wing 
Blue wing 



54 
52 



917 
467 



8585 
3804 



1013 

5871 



3 Results and Discussion 

3.1 Observed Changes in the Intensity 

The results of the conditional averaging are summarised in Figure [21 which shows the 
average intensity of the red wing measured for 54 X-class flares (in black), and for 917 
M-class flares (in gray). Time zero [t = 0] corresponds to the peak intensity of the SXR 
flux. For each flare occurring at time tf, we isolated a sequence running from tf — 90 
minutes to tf + 90 minutes and averaged it. 
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Figure 2: Conditionally averaged intensity of the red wing, measured for X-class flares 
(dark thin line) and M-class flares (grey line) 

Figure [2] reveals a small but significant increase in the intensity, which cannot be 
observed in any individual event because it is masked by fluctuations. The relative 
increase in intensity is about 0.03% for X-class flares and 0.008% for M-class flares. 
Similar results are obtained for the blue wing. We conclude that flares do have an 
influence on the sodium line and thus on the lower atmosphere. C-class flares also have 
a discernible effect on the intensity of the sodium line but their influence is too weak to 
be properly studied and so we shall not consider them further. 

GOLF measures the intensity of both wings of the line, integrated over the full solar 
disk. As the solar flare occurs over a small area on the Sun (less than 1% of the total 



surface) , the resulting intensity is the sum over the intensity originating from in the quiet 
Sun, which is by far the most important contribution, and the intensity originating from 
the flaring area. This explains why the increase in intensity is relatively weak. 

3.2 What Causes the Intensity to Change? 

Several mechanisms can be invoked to explain the observed intensity change during 
flares. The main ones are: 

• The temperature minimum region and the upper photosphere are heated, which 
results in an excess continuum emission. This emission can originate from Paschen 
and/or H~ contin uum. This effect has been invoked to explain WL emissions 
during solar flares (jXu et all 120061 : iHudson. Wolfson. and Metcali 12006 ) . We then 



expect the entire intensity profile of the sodium line to increase and so both wings 
should vary by the same amount. 

• Sodium is heated and excited or ionised. This causes a depletion of the absorbing 
sodium and hence an increase in the intensity of the absorption line. The resulting 
signature is the same as before. 

Hydrogen lines are know n to change som etimes from absorption to emission during 
very large X-class flares ( Svestkal . Il976l ). These lines, however, are mostly gener- 



ated in the upper chromosphere where more energy is released than in the upper 
photosphere. As the Na D lines are generated in the upper photosphere, we expect 
them to remain in absorption for all the flares that we consider. 

• Sodium is heated and some material is rising. This causes both a Doppler shift 
and an intensity change. The upward flow of sodium shifts the line towards the 
blue and so the intensity decreases in the blue wing with respect to the red one. 

Other effects may contribute to this. The recombination of chromospheric sodium 
ions, for example, which is enhanced by the flare-associated density increase, may shift 
the sodium line to the blue and hence lead to a decrease in the blue wing. 

If the GOLF sodium line had been measured simultaneou sly at different wavelengt hs 
(as will be the case with the proposed GOLF-NG instrument ( Turck-Chieze et alU200w ). 



then a partial separation between these different mechanisms will be possible. Our ob- 
servations, which are based on the measurement of one single wavelength at a time, a 
priori totally exclude such a separation. The important point, however, is that even 
though we do not measure the red and blue wings simultaneously, we do observe their 
average temporal evolution under similar conditions. For that reason, the two wings 
are directly comparable in a statistical sense even though individual events are never 
recorded simultaneously. A few assumptions, however, need to be made before we can 
pursue with such a comparison. 

First, our events occur at different times, so we have to assume that the impact of 
flares on the photosphere is not solar-cycle dependent. We can indeed reasonably assume 



that the speed of plasma motions should be about the same, regardless of the phase with 
the solar cycle. 

Secondly, for our samples to be statistically comparable, the number of events and the 
average SXR flare intensity should be about the same in both categories. We checked 
this and the sizes of the samples, which are given in Table [H are indeed sufficiently close. 
In particular, the noise level, which is proportional to l/yiV (where N is the sample 
size) should be approximately equal in both wings, for a given class of flares. 

Since GOLF observes the full Sun, the sodium line integrates contributions coming 
both from the quiet Sun and from active regions. Because of the statistical procedure, 
however, the former are averaged out whereas the latter add up coherently. Obviously, 
the observed variation in the integrated intensity is considerably smaller than that ob- 
served in the vicinity of the flare. 

Another point is the location of the solar flare on the Sun. Because of projection 
effects, flares occurring at the center of the disk and at the limb have different Doppler 
shifts along the line of sight. The shift associated with radial motion is largest for flares 
located close to the center of the disk. To avoid such geometric effects, we restricted our 
sample to flares whose heliocentric angle is less than 30 degrees. 

3.3 Interpretation of the Intensity Increase 

Let us start with X-class flares, whose signature on the sodium line is most evident. 

Figure [3] shows the conditionally averaged intensities of both wings, together with 
their difference {e.g Led - Ihuo). The two samples have similar sizes (see Table [T]) and 
similar average X-rays fluxes, so a statistical comparison is meaningful. The standard 
deviations for both wings are very close. Moreover, the slow degradation of GOLF im- 
plies that instrumental errors should be about the same in both samples. 

Interestingly, the blue and red wings evolve differently in time. The intensity of the 
red wing is clearly higher at the flare maximum as compared to the blue wing. After the 
maximum, the reverse occurs, and the intensity of the red wing falls off more rapidly. 
Our interpretation is sketched in Figure HI The line represented in this sketch can be 
either Na Di or Na D2, as we assume that they behave similarly. 

First, the intensity of the whole line profile rises, due to the excitation and heating 
of the sodium. The intensity of both wings indeed increases during a flare. For that 
reason. Figures 4a, 4b, and 4c present a higher profile than for the quiet Sun. Moreover, 
the intensity difference at the peak of the fiare suggests that the sodium profile is asym- 
metric. A blueshift occurs, which means that sodium is rising, see Figure 4a. After the 
peak of the flare, a deficit in intensity occurs in the red wing, which implies a downward 
motion of the plasma. A redshift may also occur, as represented in Figure 4b, meaning 
that the plasma is falling. Finally, the excitation diminishes, see Figure 4c. Note that 
the lack of observation of the core of the line prevents us from concluding about global 
changes in the line profile. The pertinent result, however, is the difference between both 
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Figure 3: Top panel: conditionally averaged intensity of the red (in gray) and the blue 
(in black) wing for the X-class flares. Curves have been smoothed over two 
minutes. Bottom panel: intensity difference between the two wings (Led - Ibiue) 

wings, which can be unambiguously attributed to a Doppler shift. 

The intensity of both wings comes back to pre-flare values in approximatively 100 
minutes, as shown in Figure [3l Not surprisingly, this duration is close to that observed 
in the SXR flux. Figure [His meant to be schematic only. We do not know, for example, 
what is really happening in the core of the line. 

Similar conclusions apply to M-class flares (Figure [5]), with a lower signal-to- noise 
ratio and a shorter recovery time. The difference between the red and the blue wings 
(bottom panel) is not properly resolved during the impulsive phase of the solar flare, 
while the redshift is clearly visible. These results provide direct evidence for the impact 
of flares on the photosphere and not just on the chromosphere. 



Photospheric sodium is affected by several effects during a solar flare. All of these 
effects (excitation, heating, motion) can change the profile of the Na D lines. The speed 
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Figure 4: Schematic representation of the possible causes for the intensity changes ob- 
served in Figure 3. The triangles represent the absorption profile at rest, (a) 
blueshift: Upward motion to chromosphere, (b) redshift: Plasma is falling, (c) 
Reduced excitation of sodium: whole profile rises down, always with a redshift. 

of ascent could be determined using the Doppler shift caused by the plasma motion. 
These effects however are simultaneous, which means that the change of intensity is 
not just due to the Doppler shift. Furthermore, we sample the line at two wavelengths 
only, which is not enough to distinguish, if they were simultaneous the contribution from 
heating from that of the Doppler shift. 



Our interpretation of the plasma motion, with first an upflow and then a downflow of 
the sodium during the flare, in based on the hypothesis that the Na D lines remain in 
absorption. We already discarded flares whose peak energy exceeds 5 x 10~^ Watts m~^ 
in order to exclude the most energetic events. Since the behaviour of the lines seems to 
be the same for X-class and M-class flares, one can reasonably assume that the Na D 
lines always remain in absorption for the events under consideration. 

These results are a priori in contradiction with the observations of downflow in the 



Ha line (JFisheiJ . Il987al ). Indeed, in the classical picture of chromospheric evaporation, 
the energy deposition layer is located in the middle of the chromosphere, and so one 
would expect photospheric lines such as Na D to behave similarly to the Ha line. A 
reconciliation between these two pictures is difficult unless: 

• The sodium affected by solar flares is not the photospheric sodium normally mon- 
itored by GOLF for a quiet Sun, but the recombined chromospheric sodium that 
emits during solar flares. The lines observed by GOLF should then behave like a 
chromospheric one. The disagreement with the Ha lines then is not so surprising 

• The affected sodium originates from the photosphere. Our results then suggest that 
at least part of the energy can reach lower altitudes than expected in the classical 
p icture. The upper photosphere c a.n indeed b e heated by several mechanisms 
( (JAboudarham and Henouxl . Il986l ). (jDingi . 120071 )) leading to white light emission 
from both the chromosphere and the photosphere. Recent studies strongly suggest 
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Figure 5: Top panel: conditionally averaged intensity of the red (in gray) and the blue 
(in black) wing for the M-class flares. Curves have been smoothed over two 
minutes. Bottom panel: intensity difference between the two wings. 



that these white -light flares are more common than expected (jFletcher et aLl . 120071 : 
Jess et g/.l . I2OO8I ) . As we do indeed analyze the average over several flares, and not 



a particular event as for the Ha line, white-light flares could have a predominant 
weight in our statistics. 

These results and their interpretations are based on only two points in the spectral 
profile, so the conclusion on the sodium dynamics in the lower solar atmosphere needs 
to be confirmed by further studies. Our results nevertheless strongly suggest that the 
upper photosphere is affected by solar flares. 

4 Conclusions 

The aim of this study was to investigate whether solar flares impact the Na Di and Na 
D2 Fraunhofer lines, as measured by GOLF. A statistical analysis based on conditional 
averaging reveals intensity increases during X and M-class flares, which demonstrates 
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that the Na D Unes are affected. Although GOLF observes only one wing at a time, the 
average evolution of the red and blue wings can be compared in a statistical sense. The 
observed intensity differences suggest the existence of plasma motions. The blueshift that 
coincides with the peak of the flare can been interpreted as a upward plasma motion, 
which is most likely associated with photospheric Na rising towards the chromosphere. 
The peak of the flare is followed by a redshift, which corresponds to a downward flow. 
Since the Na D line is mostly generated in the upper photosphere, our study thereby 
provides strong support for the influence of flares on the lower solar atmosphere. 

This interpretation is a priori in contradiction with the classical picture of chromo- 
spheric evaporation. It can however be explained if we assume that white-light flares, 
for which the lower solar atmosphere needs to be heated, have a preponderant effect 
in our statistics. White-light flares could indeed be more common than expected, and 
our study supports this idea. However, further quantitative studies are necessary before 
we can conclude that solar flares do really affect the lower solar atmosphere. Plasma 
motions could be more accurately determined with more than two measurements on the 
profile of Na D lines. The GOLF-NG has been proposed for that purpose. Another 
important conclusion is that future helioseismic instruments should take into account 
the impact of solar flares in order to improve the measurement of solar oscillations. 
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